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Abstract 

(n) 

Let Pj be the probability that a Hadamard quantum walk, started at site j on the 

integer lattice {0, . . . , n}, is absorbed at 0. We give an explicit formula for pi 1 ■ Our 
formula proves a conjecture of John Watrous, concerning an empirically observed 
linear fractional recurrence relation for the numbers p± . 
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1 Introduction 



Consider a Hadamard quantum walk on the sites 0,1,..., n, as denned in 
[I]. The boundary sites, and n, are absorbing, so any walk is certain to be 

(n) 

absorbed. Let pj denote the probability that the walk ends at location 0. 

The main result of this paper is an explicit formula for this probability. Let 
A = 2 + \/2, and B = 2 — \/2. We will prove that when n > 1 and 1 < j 1 < n, 
we have 

(B) _ V2 (A»-i - W 1 + g^j 

Pj 4 A n ~ l + B n ~ l ' 
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In particular, when j = 1, 



72 A"- 1 -j?"' 1 

Pi = 1 r, (1-2) 

2 A" -1 + B n ~ l 

and this satisfies the recurrence relation 

,„, i + 2 P 'r " 

This recurrence relation, conjectured by Watrous from numerical data, appears 
in [lj. It is a consequence of (I1.2D that 

a result that was proved with some effort in p]. 

The absorption probability values are interrelated in many interesting ways. 
In the two-dimensional table of pf \ there is a linear fractional recurrence re- 
lation, similar to (11.31) . for each column and each diagonal. This "numerology" 
was first observed empirically, and then led us to conjecture (11.11) . There is 
also a linear recurrence relation common to all rows: if 1 < j < n — 3, we have 

Two other relations, which we discuss later, have combinatorial interpreta- 
tions. 



It is interesting to consider the implications of (11. ip for starting sites in the 
"interior" of the lattice. For j = an, a fixed and n — > oo, Pj will be close to 

the limit ^ = 0.35355... . Thus, from the interior region, the probabilities for 
absorption at the left and right are almost constant, approximately 35% and 
65%. On the other hand, for the classical random walk, with equal probability 
of moving left and right, the probability that the walk, starting from j, is 
absorbed at the left, decreases linearly with j, from 1 at the left barrier to at 
the right barrier. (One proof of this appears in |3jj.) This gives another example 
of the idea that quantum walks "spread out" more evenly than classical walks 
do. 



2 Some Generating Functions. 



Our work will be based on the path count generating function approach which 
was employed in p] and [2]. 



2 



Recall that for our walk, each site j has two states, corresponding to the walker 
facing left and facing right. The one-step evolution matrix is a Hadamard 
transformation, i.e. starting from state \n, R) the particle can go next to 

\n + 1,R) with amplitude 1 / \Pl 
\n — 1, L) w. a. 1/V2 

and from \n, L) it can go to 

\n+ 1,R) w.a. 1/V2 
\n — l,L) w.a. —1/V2 

The initial state of the particle is \j, R), with < j < n. 

We define the path count generating function to be 

/j n) (*) = £ fe<WU m (2.1) 

m>\ \ P / 

where the m-th sum is over paths P of length m that are absorbed at the left 
(0) state, and the sign, cr(P), is (—1) ^ °^ blocks . In computing the sign, 
overlaps count, for example, the sign of LLL is +1. 

The probability of absorption at is given by the formula 

Pf = ^/, JifWlVMz, (2.2) 
J 2m J\z\=i/V2 J 

In pQ it was shown that 

f? = z (2.3) 

and for n > 1, 

/W = Z ±—Wi. (2.4) 

We now complement this with another recurrence relation that allows (12. ip to 
be computed for j > 1. Observe that any path from j that is absorbed at 
must go through 1. Hence any such absorbed path breaks up into: a) a path 
from j to 1 reaching 1 only once; and b) a path from 1 to 0. Part a) is of the 
same shape as a path from j — 1 to on a lattice with absorption at n — 1, 
and part b) is just a path from 1 to 0. The path of part b), if it immediately 
moves left, must be preceded by an L move, so we must correct the sign for 
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this case (and this case only). This gives 

ff = ffr? (A (n) - 2z) , 2<j<n. (2.5) 



Here are few of these functions. 



f 0) z{2z* - 1) (3) _ 



(z-i)(z+iy « (z-i)(z+i) 

lU z{Az 4 - 3z 2 + 1) (4) _ z 2 {2z 2 - 1) (4 ) = ^ 3 

Jl 22 4 -2z 2 + l ' /2 2z 4 -2z 2 + l' /3 2z 4 -2z 2 + l' 
AS) = z{2z 2 - \){Az 4 - 2z 2 + I) (5) 

1 /'o..?! i ..9 .. 1\/n...1 -.9 .. i 



z(2z 2 -l)(4z 4 -2,z 2 + i) (5) ^ 2 (4^ 4 -3^ 2 + l) 

(2z 3 + 2 2 - ^ - l)(2^ 3 - ^ 2 - ^ + l) ' * 2 ~ (2z 3 + z 2 -z- l)(2z 3 - z 2 - z + l) ' 



f (5) = ^ 2 ~ I) (B) _ 



(2z 3 + ^ 2 - z - 1)(2^ 3 - z 2 - z + 1) ' J4 (2^ 3 + z 2 - z - l)(2z 3 - z 2 - z + 1) ' 

We note that for each n, the /• have the same denominator. This can be 
proved as follows. We first observe that the power series for fj begins with ±z J , 
so that fj = z^u, with u(0) ^ 0, oo. In particular, we have fi = za(z)/b(z), 
where a and b are polynomials. Next, we group paths according to the location 
of the first R move and find the relation 



fj = z fj+l + + I - )' ' : 



k=2 



valid for 1 < j < n — 1. From this it follows that is a polynomial, using 
induction on j. 

It is interesting that the path sign is essentially the same as the Rudin-Shapiro 
coefficient. This coefficient a n is determined by the parity of the number of 
"11" blocks in the binary notation of the positive integer n. The Rudin-Shapiro 
coefficient has many applications, including the solution of extremal problems 
in classical Fourier analysis. For a survey of this topic, see jl]. 



3 Some Combinatorial Results. 



It is interesting to see how much can be determined by purely combinato- 
rial arguments, without relying on integration. We begin with two interesting 
relations. 

Theorem 3.1 We have 

„(«) i „(«) _ i 

Pi + Pn-l — l - 

That is, in any row, the outer entries sum to 1. 
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Proof. Let P be a path that starts from 1 and is absorbed at n. Its complement 
P, obtained by interchanging L and R, is a path that starts at n — 1 and is 
absorbed at 0. We note that P must begin and end with R moves. Let P 
contain i L moves, r R moves, and have k occurrences of RL. Then the sign 
of P is ( — iy~ k . However, 

k= # of RL in P = # of LR in P , 

since P begins and ends with L moves. Therefore we have 

a{P) = {-iy~\ a{P) = (-iy- k , 

which implies 

a(P)a(P) = (-lf +r = (-l) r - £ = (-I) 71 ' 1 . 

This tells us that if we complement all paths, the probability is unchanged, 
since it is a sum of squares of quantities (signed path counts) that individually 
change only by a sign. So 

1 — Pi = Pr[ a walk from 1 is absorbed at n } 
= Pr[ a walk from n — 1 is absorbed at ] = p„_i- 

□ 

If we let qj 1 ^ = 1 — pj be the probability that the walk reaches site n, then 
Theorem 13.11 looks like an "obvious" symmetry relation that should also hold 
for j > 1. However, this is not so. (See Table 1 at the end of this paper.) 

Theorem 3.2 We have 

2p£ n) =;4 n) + l. 

In words, doubling the first number in any row is the same as increasing the 
second by 1. 

Proof. Observe that 

fi n \z) = z + zff\z). 
Now take the Hadamard square of both sides, and evaluate at 1/2. The two 
terms in the right side do not interfere because f^i?) is a multiple of z 2 . □ 

For n — 2, 3, the absorption probabilities can be obtained by combinatorial 
reasoning, without any integration. First, for n = 2, we have 

pf ] = 1/2, 

since there are only two possible paths, each absorbed after one step. To 
compute p^\ observe that for t = 1, 3, 5, . . ., there is precisely one path that 
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reaches site after t steps. Therefore, the signs are irrelevant, and we have 



(3) 
Pi 



111 2 
2 + 8 + 32 + "'~3' 



Using Theorem 13.11 we find 



(3) 
P2 



1 

3' 



4 Proof of the Explicit Formula. 



In this section we prove that (jl.ip holds. What we would like to do is integrate 
(12. 2p by residues, and expose the dependence of p^ on n. With the original 
integral, this is probably impossible, since we do not know where the poles 
inside the circle of integration actually are. However, we can express p^ using 
the integral of a new rational function, with the same mysterious poles inside, 
but with only one pole outside. Since the sum of the residues of any rational 
function vanishes, we can just as well evaluate the residue outside the circle, 
and this leads to a formula for Pj. 

Let us begin with another formula for the path count generating function. Let 

a, j3 be the two roots of 

T 2 — (1 — 2z 2 )T - z 2 = 0. 

These will always be used symmetrically so we do not care which is which. 
Explicitly, 

n l-2z 2 ± VI + 4z 4 

,H 2 
Also there is a recurrence relation 



a: 



■ k - f l -2z 2 )a k ~ l + z 2 a k ~ 2 (4.1] 



and similarly for (3. 

Lemma 4.1 If 1 < j < n, then 

fin) = z j_ 



a n-j _ pn- 



a n - (3 n + 2z 2 (a n ~ 1 - (3 r - 



-IV 



Proof. This can be proved by induction. First, let j ' = 1 and increase n, using 
(12. 4p . Then, for each n in turn, let j = 2, . . . , n — 1, and use (I2.5p . □ 

It will be convenient to allow let j = n, and define f^ = 0, so that p^ = 0. 
This is consistent with the above formula, as well as (12.51) . 
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As a function of z, fj{n) is odd or even according as j is. We now write 



f\ n ^ = z^g\ n \ and t = z 2 . Then, 



gf\t) - — 



T n 2tv n _i 
where 

a k -(3 k 
a — f3 

Since r k+2 — (a + (3)r k+ i + a/3r k = 0, we have r = 0, r\ = 1, and 

r k+2 - (1 - 2t)r k+1 -tr k = (4.2) 
Consequently, the r k are polynomials in t. 

On the circle \z\ = l/y/2, we have z = l/(2z), and t = l/(4£). Making this 
substitution into g~j, clearing denominators, and observing that a [3 = —t, we 
get 

2 

|>>| 2 = (_)'W- 



(r n + 2tr„_i)(r n - r n _ x ) 
Since = z^g^ \ we get from this 



p W = izll! / dt (4.3) 

27r? J\t\=i/2 (r« + 2tr n _i)(r n - r n _i) 



The next task is to study the poles of the integrand. 

Lemma 4.2 T/ie zeroes of r n+1 — r n are inside the circle \t\ = 1/2, and the 
zeroes of r n+ \ + 2tr n are outside it. 

Proof. The rational function is analytic for \z\ < l/y/2 (this is implicit 
in the proof of Lemma 17 of [1]), and this holds for f^ as well, since it has 

the same denominator as f[ n \ The result is a consequence of this and the 
computations used to derive (14. 3p . □ 

As a consequence of this lemma, we can choose an e > so that pf^ is given 
as in (14.31) . but with the contour of integration now \t\ = 1/2 — e. 

Lemma 4.3 Let n > 2. If 1 < j < n, H 5 = t j - l (l + 2i)r n _ 3 - + (-l) j (rj - 
r j-i)( r n + 2ir n _i) is a polynomial in t, divisible by r n — r n _i. 

Proof. The recurrence relation on rj implies that -f/j+2 = (1 — 2t)Hj + i + Hj, so 
it suffices to check the statement for j = 1 and j = 2. It is easy to check that 
H\ = r n _i — r n , and the recurrence relation on r n implies that H2 = r n — r n _\. 
□ 
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The above proof has the consequence that Hj/ (r„ — r n _i) depends on j only. 

Theorem 4.4 Letn>2 and 1 < j < n. Then / TO]) holds. 

Proof. Using Lemma [4.31 to rewrite t J_1 r n _j, we have for 1 < j < n 

M \_ f r n-A r i - rj-i)dt 1 f Rj(t)dt 



2m y|*[=i-e (1 + 2t)(r„ - r n _i) 2ttz ./|t|=A-e (r n + 2tr„_i)(l + 2t) ' 

for some polynomial Rj. The second integral is zero, because all the poles of 
the rational function are outside the integration contour. In the first integral, 
the integrand 

r n~j( r j ~ r j-l) 



;i + 2t)(r n -r n _ i ; 



as a function on the Riemann sphere, has a unique singularity outside the 
contour, which is a pole of order 1 at t — — |. Indeed, the degree of r\ is k — 1 
so the degree of the denominator is two plus the degree of the numerator, 
which assures that there is no pole at infinity. 

Since the residues of a rational function sum to zero, we get 

(n) = R i rn-fa-r^) = l r -(rj-r Q / IN 
^ 4 — 3(l + 20(r n -r n _ 1 ) 2 (r n - r n _0 V 2/ y } 



2 



To derive the explicit formula, observe that a ( — 1/2) = and /3(— 1/2) 
^, so (with A, ^ = 2± y/2) 

1 A fe — 1 A k ~ 1 — 

r fc(~2'' = 2 fc+1 /2 ' ( r fc - r fe-i)(--) = • 



Then, substitute these values into (I4.4p and simplify. □ 

It is natural to extend our notation so that = 1 for n > 2. The explicit 
formula does not work there, but the above proof indicates a reason for this. 
We could use (14.21) to extend rj to j = —1, but then r_i would not be a 
polynomial, invalidating our arguments. 



Using the recurrence relation for r&, we can prove that when n > 1, the 
polynomial r n — r n _i has distinct roots. However, the computations for this 
are not very enlightening, so we leave verification of this to the reader. 
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5 Remarks. 



Initially, we had arrived at the formula (II. II) after numerical calculations that 
were done by a different method, which we believe to be of independent in- 
terest. In this section, we discuss how these calculations were done. The idea 
is to combine numerical approximation of the residues with a bound on the 
denomininator. 



Suppose we want a value for 

in which / G Q(z), and a G Q. Let a be large enough that all the poles of / 
are inside the circle \z\ = a. On this circle, z = a 2 /z, and if we use this to get 
an expression for /, we can bring (15. ip into the form 

mi 1 f b(z)dz 



p = 



m 2 2iri J\ z \=a c(z)d(z) ' 

in which mi,m 2 G Z, b, c, d G Z[t>], and c,d are monic. The zeroes of c and 
d are algebraic integers, and we choose notation so that those of c and d are 
outside and inside the circle, respectively. 

In the cases of interest to us, c and d had distinct roots, so let us make 
this simplifying assumption. Then, evaluating (2ni)~ 1 J b(cd)~ 1 dz by residues 
produces 

In this expression, £ and £' range over zeroes of d, and rj ranges over zeroes of 
c. The denominators are algebraic integers, and 

TT(£ ~ ^) I resu hant(c, d) := R; 
v 

II (£ — CO I discriminant (d) := D. 

Since 6 has integral coefficients, we conclude that ( 15. ip is an integral multiple 
of A -1 , where 

A = m 2J R£>. 

In particular (by Galois theory), p is a rational number. 

To obtain (15.11) exactly, then, it will suffice to compute the integer 5, and 
then evaluate (15. 2p using numerical approximations to the zeroes, with enough 
accuracy to determine (15. ip to the nearest integer. 
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Using this method, we were able to compute absorption probabilities exactly 
up to n = 20 in a couple of minutes on a workstation. Straight numerical 
integration would have been much slower, and would not have given us exact 
results. 

We end this paper with a short table of the Pj. The numerators of p\ are 
Sequence A084068 in [5]. 

TABLE 1. Absorption Probabilities Pj ■ 





3 = 1 


2 


3 


4 


5 6 


2 


1/2 










3 


2/3 


1/3 








4 


7/10 


4/10 


3/10 






5 


12/17 


7/17 


6/17 


5/17 




6 


41/58 


24/58 


21/58 


20/58 


17/58 


7 


70/99 


41/99 


36/99 


35/99 


34/99 29/99 



8 239/338 140/338 123/338 120/338 119/338 116/338 99/338 

9 408/577 239/577 210/577 205/577 204/577 203/577 198/577 169/577 
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